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For a semiconductor quantum dot strongly coupled to a microcavity, we theoretically investigate 
phonon-assisted transitions from the exciton to a cavity photon, where the energy mismatch is 
compensated by phonon emission or absorption. By means of a Schrieffer-Wolff transformation we 
derive an effective Hamiltonian, which describes the combined effect of exciton- cavity and exciton- 
phonon coupling, and compute the scattering rates within a Fermi golden rule approach. The 
results of this approach are compared with those of a recently reported description scheme based 
on the independent Boson model [U. Hohenester et al., Phys. Rev. B 80, 20131 1(R) (2009)], and a 
numerical density matrix approach. All description schemes are shown to give very similar results. 
We present results for the spontaneous emission lifetime of a quantum dot initially populated with 
a single exciton or biexciton, and for the spectral properties of an optically driven dot-cavity system 
operating in the strong coupling regime. Our results demonstrate that phonon-assisted feeding plays 
a dominant role for strongly coupled dot-cavity systems, when the detuning is of the order of a few 
millielectron volts. 
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I. INTRODUCTION 

Cavity quantum electrodynamics (QED) investigates 
hybrid systems consisting of a quantum emitter coupled 
to a quantized cavity mode. It has been first imple- 
mented with atoms in high- finesse optical resonator s.^^^^ 
Semiconductor systems allow for the design of monolithic 
structures."^ A particularly promising approach is based 
on semiconductor quantum dots embeddedjn microcav- 
ities, such as micropillars ,^*^*^ ' ^ * microdisks,^ or photonic 
crystal nanocavities.^^ With theses systems many cav- 
ity QED phenomena have been observed, including Pur- 
cell enhancement,^ photon bunching and antibunch- 
ing,^ single and entangled photon sources,^ and 
strong light matter coupling .I^ESEIl Possible applications 
range from ultralow threshold nanolasers,^over efficient 
single and entangled photon sources, to various applica- 
tion in quantum information science.'^' 

In contrast to atoms, semiconductor quantum dots, 
sometimes referred to as artificial atoms^ are in intimate 
contact with their solid state environment. This leads 
to enhanced scattering and dephasing contributions. In- 
deed, cavity-QED experiments with semiconductor quan- 
tum dots have revealed a significant feeding channel for 
the cavity mode,*^"^^*^ even for largely detuned dot- 
cavity systems, which has been attributed to the influ- 
ence of the environment. However, the microscopic ori- 
gin of such feeding has been discussed controversially in 
the literature. Shake-up processes in charged quantum 
dots^^ and quasi-continuum excitonic transitions^^ have 
been suggested as possible processes for populating the 
cavity mode, even in presence of large detuning of several 
tens of meV. For smaller detunings, of the order of a few 
meV, phonon processes are expected to play an impor- 
ta nt rol e. It is now well established that phonon dephas- 
in^^^ and relaxatio n^^ ' ^^"^^ govern the coherent optical 



response of quantum dots. As for cavity-QED, theoreti- 
cal work has suggested that such pure phonon dephasing 
might be responsible for cavity fee ding at detun ings of 
the order of the polariton linewidths )^^'^^ ' ^'^'^ -^^ Experi- 
mental work has confirmed th at for such small detunings 
dephasing plays a crucial role.^^I^ 

In a recent paper, it has been demonstrated theoreti- 
cally and experimentally that, besides phonon dephasing, 
there exists another phonon-mediated feeding channel for 
the cavity that works efficiently over a broad detuning 
range of several meV. Similar results have been indepen- 
dently obtained by other author s. l^^^ The cavity feeding 
relies on a combined effect of exciton-cavity and exciton- 
phonon couplings, and accounts for a process where the 
exciton decays into a cavity photon and the energy mis- 
match is compensated by the emission or absorption of a 
phonon. 

In this paper, we present the detailed analysis of the 
theoretical framework given in Ref . ^32 , and provide a sim- 
ple description scheme for phonon-assisted cavity feed- 
ing. More specifically, we bring, by means of a uni- 
tary Schrieffer-Wolff transformation, the electron-phonon 
hamiltonian to a form that precisely describes the com- 
bined effect of exciton-cavity and exciton-phonon cou- 
plings. We compare our results with those of the previ- 
ously used independent Boson model and of a numerical 
density matrix approach, and find almost identical re- 
sults for all approaches. We also present results for the 
decay of excitons and biexcitons in presence a phonon- 
assisted cavity feeding, and analyze the impact of such 
feeding on the spectral properties of a driven dot-cavity 
system. 

We have organized our paper as follows. In Sec. [ll| we 
present the theoretical description scheme for phonon- 
assisted cavity feeding in semiconductor quantum dots. 
We start with the Schrieffer-Wolff transformation, where 
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the cavity and phonon couplings are treated as pertur- 
bations, which ahows us to compute the feeding rates 
by means of a simple Fermi-golden-rule approach. Due 
to the simple structure of the independent Boson hamil- 
tonian, one can also describe the phonon dynamics ex- 
actly and only treat the exciton-cavity coupling as a per- 
turbation. We finally introduce a density matrix ap- 
proach, that allows for simulations without employing 
any approximation. Some details of this approach are 
presented in the appendix. In Sec. [ill] we present re- 
sults for these different description schemes. We discuss 
the feeding rates and the resulting spontaneous emission 
lifetimes for different detunings, temperatures, and cav- 
ity quality factors. We also investigate the spectra for a 
driven dot-cavity system, and analyze the detuning de- 
pendence of the polariton intensities. Our results demon- 
strate that phonon-assisted feeding plays a dominant role 
for strongly coupled dot-cavity systems when the detun- 
ing is of the order of a few meV. Finally, in Sec. IV we 
summarize and draw some conclusions. 



II. THEORY 



A. Model 



Hamiltonian 



In the following we consider a quantum dot embedded 
in a microcavity, which interacts with phonons.*^^-^ We 
will show that, due to the coupling of the quantum dot 
state with the phonons, excitonic states can decay into 
cavity photons even for relatively large exciton-cavity de- 
tunings of a few meV. We start with a description of the 
excitonic states in terms of a generic two level system, 
where e denotes the excited and g the ground state. Al- 
though we shall be primarily concerned with the exciton 
to groundstate decay, our description scheme also applies 
to the biexciton to exciton decay. With the quantum dot 
operators agg = \g){g\ and age = |e)(e|, and a and bk the 
bosonic annihilation operators for the cavity mode and 
the kth. phonon mode (with energies cjcav and cj/^), the 
Hamiltonian for the uncoupled system reads 



Ho = {EgCFgg ^ E ^(7 ^ ^ C3M ^ ^ k ^l^fc ' 

k 



Here Eg and E^ are the energies of the ground and ex- 
cited quantum dot state. We set = 1 throughout this 
paper. The exciton-cavity coupling in the rotating- wave 
approximatiori^^^ is of the form 



Hg = g {(Jega + (Jg 



(2) 



where cFge = \g){e\ and a eg = \e){g\ are the lowering and 
raising operators for the quantum dot, and g is the cou- 
pling constant, which depends on the exciton dipole mo- 
ment and the field distribution of the cavity mode.'^^'^ 
Through Hg exciton and cavity become coupled, and new 



polariton modes are forme d which have partial exciton 
and partial cavity char act er.'^'^^'^ In the strong coupling 
regime, the coupling is stronger than the quantum dot 
and cavity losses, and excitation is coherently transferred 
between the dot and the cavity. In contrast, in the weak 
coupling regime losses dominate over the coherent cou- 
pling, and the excited state decays mono-exponentially. 
Quite generally, the calculation of the phonon-assisted 
scattering rates, to be discussed below, works in both 
the strong and weak coupling regime, and we thus will 
not be specific about this point until Sec. |II P] where we 
introduce our master-equation approach. 

The coupling between the quantum dot states and 
the phonon s is desc ribed by the generic electron-phonon 
HamiltoniarPE^ 



H, 



ep 



^99-^99 



' ^ee-Bee i 



where we have introduced the hath operatorJ^ 



Bj 



k 



(3) 



(4) 



Here, \ij^k denotes the electron phonon matrix elements. 
In case of the exciton to groundstate decay Bgg is zero, 
since only the exciton interacts with the lattice degrees 
of freedom, whereas for the biexciton to exciton decay 
both Bgg and Bee are non-zero. The diagonal bath oper- 
ators Bgg and Bee describe a coupling where the lattice 
becomes distorted, but no transitions between quantum 
dot states are induced. The corresponding mode l is us u- 
ally denoted as the independent Boson model?^ 35 |39|4Q |Jj^ 
Eq. Q we have defined the bath operators in a slightly 
more general way, such that also transitions between dif- 
ferent quantum dot states can be described, for reasons 
that will become clear further below. 



2. Scattering rate 

To compute the phonon-assisted scattering rates be- 
tween the quantum dot and the cavity, we employ a mas- 
ter equation approach. Let 



H = Ho^V 



(5) 



be the Hamiltonian composed of a part Hq, which can be 
treated exactly, and a perturbation part V. In the most 
simple approach, is associated with the Hamiltonian 
of Eq. ([T]), and V with the exciton-cavity and exciton- 
phonon couplings of Eqs. Q and ([3|. However, as we 
will show in Sec. |II C[ due to the simple structure of the 
independent Boson hamiltonian of Eq. ([3| we can even 
describe the full dynamics of Hq + Hq^ exactly, and only 
consider Hg as a perturbation. 

The time evolution of the composite quantum-dot cav- 
ity system, described by the density matrix p, and the 
phonon part pph, can be approximately described by the 
master equation in Born approximation^^^I^ 

dpjt) 
dt 



^ trph ( [V{t), [V{t), p{t) (S> /Oph] ])dT.^ 



(6) 
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where V{t) is the perturbation in the interaction rep- 
resentation of In Eq. (|6| we have assumed a fac- 
torizable density matrix^ at the initial time to. In the 
fohowing we consider transitions between the initial and 
final states i and / of the quantum-dot cavity system, 
which, e.g., for the phonon-assisted decay read |e;0) and 
1^; 1). Here, the first entry corresponds to the quantum 
dot state and the second one to the number of photons 
in the cavity. Upon explicit evaluation of the double- 
commutator in Eq. ^ we encounter contributions where 
both V operators act from either the left or right hand 
side on p, which can be associated with generalized out- 
scatterings and dephasing.l^IEII xhe other contributions 
can be associated with generalized in-scatterings. More 
specifically, we find for the increase of population in the 
final state 

Pff = I trph [{f\V{mpuPv^(i\V{r)\f)+c.c.] dr , 

(7) 

where c.c. denotes the complex conjugate of the preced- 
ing term. We have used the adiabatic approximation of 
letting the lower integration limit approach minus infinity 
and introducing a small damping term inside the integral 
(not shown), which allows us to perform the time inte- 
gration prior to the phonon trace. From Eq. ^ we can 
infer the scattering rate 



is achieved through the operator 



Tif = 2 Re 



r trpJpph(^inT)|/)(/inO)K)) 

J —CO 



dr (8) 



for a transition from the initial state i to the final state 
/. To arrive at our final expression, we have exploited 
the cyclic permutation of operators under the trace. The 
generalized Fermi golden rule expression of Eq. ([8| will 
serve us as the starting point for the calculation of the 
phonon-assisted scattering rates. 



B. SchriefFer— Wolff transformation 

The decay of an exciton to the cavity under off- 
resonant conditions requires both the coupling between 
states e and ^, by means of the exciton- cavity hamilto- 
nian Hg^ and the emission or absorption of a phonon, 
in order to compensate for the energy mismatch. The 
scattering rate will thus depend on both cavity and 
phonon couplings g and A^. A convenient way to ac- 
count for the combined cavity-phonon character of the 
scattering is to remove, in lowest perturbation order, the 
exciton-cavity coupling by means of a Schrieffer- Wolff 
transformation.!^^™ Let us consider first the Hamiltonian 
H = Hq + Hg composed of the free part and the exciton- 
cavity coupling. By means of a unitary transformation^^ 



H 



H 



(9) 



we can remove the exciton-cavity coupling to the lowest 
order of ^, by choosing s such that [5, Hq] -\-Hg = 0. This 



S = ^ [(Jega 



(JneCl 



(10) 



where A = Ep 



is the detuning between 



the quantum dot transition and the cavity. Then, H = 
Hq + 0{g'^) describes the uncoupled exciton-cavity sys- 
tem, and the second-order corrections can be identified as 
Lamb and Stark shift which will be neglected below. 



1. Diagonal electron-phonon coupling 

Let us investigate how the unitary transformation of 
Eq. ([9| affects the electron-phonon coupling. After some 
simple calculation we arrive at 
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{<Jega - a-gett^) [Bee " Bg 



(11) 



for the lowest order correction of Hep^ which is the ex- 
pression we are seeking for. The second term on the 
right-hand side describes a combined interaction of Hg, 
which couples the quantum dot states e and g via emis- 
sion or absorption of a cavity photon, and i^ep^ which can 
compensate for the energy mismatch A via emission or 
absorption of a phonon. In what follows, we assign the 
perturbation V to the Hamiltonian given by Eq. (11), 
and consider an interaction representation according to 
Hq. Upon insertion of the Hamiltonian of Eq. (11) into 
the scattering rate of Eq. ([8|, we obtain for the transition 
from the initial state |e; 0) to the final state \g] 1) the rate 



Le^ = 2Re 



/o 



i(A-iO+) 



^ {SB{t)SB{0)) dr 



(12) 

with 0+ being a small positive quantity and SB = 
Bee — Bgg the bath fluctuation operator. We have used 
(•) = trph(pph-)- thermal equilibrium only the phonon 
occupations {h\h^) = fik are non-zero {fik is the Bose- 
Einstein distribution at temperature T), and we can eval- 
uate the phonon correlation function analytically.^^ We 
then get 

Te^ = 27r^ |Aee,fc - \gA^ 

k 

X (nfe(5(A + ujk) + (fin + 1)(5(A - iJk) ) (13) 

for the phonon-assisted feeding rate of the cavity through 
decay of the exciton state. Eq. (13) is our final expres- 



sion. It accounts for the phonon-assisted exciton decay 
to the cavity, where the phonons are assumed to be in 
thermal equilibrium. 

If the exciton decays within a cavity populated with 
n photons, the scattering rate has to be multiplied by 
n (because of the cavity field operators acting on the 
Fock state with n photons). We can also derive along 
the same lines the scattering rate Tge for the reversed 
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process, where the photon is removed from the cavity 
and an exciton is created. The resulting expression is 
identical to Eq. (13) except that A has to be replaced by 
—A. This backscattering, however, is usually inefficient 
due to the strong losses of the cavity. 



e^be ^ = b -\- {X/cj) the phonon field becomes displaced 
by a constant value. Then, 



c^|6t + -) lb 



A 



2. 



-diagonal electron-phonon coupling 



^6^6 + A(6 + 6^j + ^ (16) 



Our above derivation can be also applied to the case 
where more than two quantum dot states are involved. 
Indeed, it has been demonstrated that in the cascade 
multi-exciton decay of an optically excited quantum 
dot the emission of cavity photons plays an important 
role.fli^ Since there is practically always an energy mis- 
match between the transition energies of the quantum 
dot states and the cavity photon, phonons are expected 
to play a decisive role also in this cascade. As a sim- 
ple model for a phonon-assisted cascade decay, we con- 
sider again two states e and g which are elect romagneti- 
cally coupled through Hg. However, for the phonon cou- 
pling we assume that the lower state g is coupled through 
phonon emission or absorption to some other state g^ viz. 



gives the free-phonon hamiltonian for the displaced 
phonon fields. Let us first assume that only the upper 
state interacts with the phonons, i.e., Hq^ = cieeBee (the 
case where phonon coupling is also present in the lower 
state, such as for the biexciton to exciton decay, will be 
addressed below). The transformation e^^^^ with 



(17) 



99' 



9' 9) -^99' 



(14) 



The additional state g' might be even optically inactive, 
corresponding, e.g., to the excitation of one hole or elec- 
tron out of the groundstate of a multi-exciton state. It 
turns out that for this off-diagonal phonon Hamiltonian 
the scattering rate can be computed in a completely sim- 
ilar fashion, and we arrive at 

k 

X \^nkS{A' + ujk) + {fin + l)S{A' - cj/e) ) , (15) 



then allows to generate the independent boson hamil- 
tonian through^^ Hq 4- i^ep ~ e'^^^^i^oe"'^^^*, apart from 
an energy renormalization — dee XI A^^ ^/cj/c of the quan- 
tum dot state, which is usually referred to as the polaron 
shift and which ca n be easi ly absorbed into the definition 
of the dot energies. '^^'^^'^ This transformation is known 
as the polaron transformation^ as it displaces, depend- 
ing on the state of the quantum dot, the phonon fields to 
the new equilibrium positions. For the exciton decay, the 
groundstate is associated with an undistorted lattice, and 
only for the exciton state the lattice becomes distorted, 
due to the formation of a polaron. Then, we obtain from 
Eq. ([8| the scattering rate 

.0 

Teg = 2g^ Re / trph 
J —00 

X (pph<^; lle^^'^l^; l)<e; 0|e^e-^^°^e-^ |e; O) ) , 

(18) 



-UJc 



is the detuning with respect where we have used e"" 



where A' = Ee - Eg 
to state g'. Equation (15) accounts for an indirect tran- 
sition, where the system decays via the optical channel 
to ^, and becomes promoted through phonon emission or 
absorption to the final state g\ As a result, multi-exciton 
states can easily decay even in absence of strict energy 
conservation for the exciton-cavity system. 



C. Independent Boson model 

Due to the particular structure of the independent bo- 
son Hamiltonian ([s]) it is possible to treat the electron- 
phonon interaction without further approximation. In 
this approach, we introduce an interaction representation 
according to Hq = Hq^Hq^ and treat the exciton-photon 
coupling Hg as the perturbation. 

The key to the success of the independent Boson 
mo del is p rovided by the phonon displacement opera- 
tor,™™which, for a single phonon mode, reads with 
s = X/uj{b — b^). Through the unitary transformation 



'99. 



cFee^^ ' We next evalu- 
ate all matrix elements in Eq. ([T8|) to transform the inte- 
grand to e"*^'^ (e^^^^e 



), where s{t) is the expression of 
Eq. \17\ in the interaction representation of Hq and the 



brackets denote the expectation value over the phonon 
density matrix. In thermal equilibrium the expectation 
value (e^^^'^^e*) can be evaluated analyticall}/ 4 41 , 45] 

C(t) = 

exp ( - ^ (^^y [{nu + l)e— + nu^^^^] ) . 

(19) 

This correlation function also determines pure phonon 
dephasing and the spectral li neshap e in optically excited 
semiconductor quantum dots.f^^I^ The phonon-assisted 
scattering rate within the independent Boson model then 
become^^ 



2g^ Re 



-iAr 



C(-r) dr. 



(20) 
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tion of Lindblad fornPEH 




FIG. 1: (Color online) Level schemes used in our simulations, 
(a) Ground and excited quantum dot states e, ^, and cavity. 
The exciton state is coupled to the cavity, with the cavity 
coupling g. The exciton and cavity decay with rates 7rad 
and hv. In case of laser pumping, we also consider transitions 
from g to e. (b) Biexciton level scheme, consisting of the dot 
groundstate g, the exciton states x and y, with orthogonal 
polarizations, and the biexciton state XX. In our simulations 
we use a basis of the bare quantum dot states, and these 
states together with up to two photons in the cavity. When 
the cavity and the dot are not perfectly aligned, both exciton 
transitions can couple to the cavity. We assume that the 
XX X and x ^ g transitions couple with strength gcosO, 
and the XX y and y ^ g transitions with strength g sin 
to the cavity. Here is an angle that determines the degree 
of mixing. Pumping brings the system from the ground state 
via the exciton states x and y to the biexciton state XX. 



One can easily show that, by expanding the integrand in 
Eq. (19) to lowest order, we precisely recover the result of 
Eq. (13) for the Schrieffer- Wolff transformation. Finally, 



when phonon interactions are present in both the excited 
and the ground state, we have to replace in Eq. ([19]) the 



phonon matrix elements simply by Xee,k — ^gg,k- This is 



due to the fact that only the difference in phonon equi- 
librium positions affects the phonon dephasing and, in 
turn, the phonon assisted scattering rates. 



D. Master equation approach 

With the phonon assisted scattering rates, obtained 
within either the Schrieffer- Wolff or independent Boson 
model approach, we can describe the dynamics of the cou- 
pled dot-cavity system by means of a master equation ap- 
proach^ In this work we will consider the level schemes 
depicted in Fig. [l] which consist of a different number 
of composite dot-cavity states. We describe the coher- 
ent dynamics by the Hamiltonian which accounts for 
the single-particle energies and the dot-cavity couplings. 
The incoherent dynamics, which accounts for radiative 
dot and cavity losses, phonon- assisted cavity feeding, and 
dot pumping in presence of an ex terna l laser pulse, is de- 
scribed by the Lindblad operator^^^EH , where /i labels 
the different scattering channels. The time evolution of 
the density matrix is then determined by a master equa- 



= -i[H, p]-j2(l i^i^^' p} - ^^p^i) ■ 



(21) 



Depending on the relative importance of the dot-cavity 
coupling g and the various scattering losses, the solutions 
of this master equation exhibit strong or weak coupling 
effects. The master equation of Eq. ( [21] ) can be rewritten 
in the form p = —i\-p, where L is the Liouvillian of the 
dot-cavity system. If we interpret p as a column vector, 
with matrix elements pij ordered in some unique fashion, 
L becomes a matrix, whose elements can be directly read 
off from Eq. (21). A convenient way to solve the master 



equation is to seek for the right and left eigenvectors X 
and X, which are defined through^^ 



LX = XA , XL = AX , 



(22) 



where A is a diagonal matrix containing the eigenvalues. 
The eigenvectors form two biorthogonal sets with XX = 
XX = 1. The master equation, subject to the initial 
condition p{0) = po, can then be easily solved through 



Pit) = (Xe-'^'Xj po, 



(23) 



with the expression in parenthesis being the time evo- 
lution operator for the density matrix. Eq. (23) is ex- 



tremely useful for computational purposes, since it allows 
for the full solution of the master equation by simply di- 
agonalizing once the L matrix. 

From the solutions of the eigenvalue problem of 



Eq. (22) we can also obtain analytic expressions for the 
spectrum of the coupled dot- cavit y syst em. We start 
from the Wiener-Khinchin theorem,!^^^ which, in case 
of continuous dot pumping, reads 



S{u;) (X Re 



POO 

/ e--(at(0)a(T))^^dr, (24) 
Jo 



where the brackets denote the expectation value at late 
times when the system has reached its steady state. 
Eq. (24) relates the spectrum to the two-time correla- 



tion function of the cavity fields, which can bej^omputed 
by means of the quantum regression theorem*^ Let A be 
the matrix of expectation values for the cavity field op- 
erator a between the composite dot-cavity states of our 
master equation approach. The quantum regression the- 
orem then allows us to express the correlation function 
as 



(at(O)a(r))^ 



tr ^Aji (Xe' 



(25) 

where we have implicitly assumed summation over re- 
peated indices, pss denotes the density matrix in the 
steady state, which can be computed as the long time 
limit of Eq. (23). Equation (25) has to be interpreted as 
follows. First, the operator at the earlier time, a^, acts 
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from the right-hand side on p^g, thereby creating a fluc- 
tuation in the system. Next, this fluctuation propagates 
in time, which is described by the time evolution oper- 
ator of Eq. (23). Finahy, we determine the expectation 



value for the operator at the later time, a, with the prop- 
agated fluctuation operator. Again, the knowledge of the 
eigenvectors and eigenvalues allows us to easily evaluate 
the two-time correlation function. Upon insertion of the 
correlation function into the Wiener-Khinchin expression 



for the spectrum, Eq. (24), we can compute the spectrum 



of the driven dot-cavity system according to 

S{io)^lmtv[AjJx{A-u;)-^x) {pssA^)^A ■ 

(26) 

In a very similar fashion we could also compute photon 
correlation functions. 



E. Density matrix approach 

We flnally present an approach suited for the numerical 
solution of the coupled dot-cavity dynamics in presence 
of phonon interactions, without employing a perturba- 
tion approach. We consider the situation where the sys- 
tem is initially in the excited dot state, and subsequently 
decays through radiative decay or coupling to the cavity 
photon, which then leaks out of the cavity. In absence of 
dot pumping, we are left with a generic two-level system 
consisting of the excited quantum dot state together with 
the dot groundstate and one photon in the cavity. We can 
describe the system in terms of an effective Hamiltonian 



eff 



99 



e9 



'9eJ 



(27) 



where 7 is the radiative decay rate of the excited quantum 
dot state and tz the leakage rate of the cavity. Note that, 
in contrast to the above analysis, here g denotes the state 
composed of the dot groundstate and the cavity photon. 
It is important to realize that any scattering, described 
by the loss terms proportional to 7 and hc^ brings the sys- 
tem to a state (dot groundstate without cavity photon) 
that is not coupled to our two-level system. 

The total Hamiltonian then consists of i^eff, together 
with the electron-phonon coupling of Eq. ([3| and the free 
phonon Hamiltonian. For simplicity, we consider the sit- 
uation where only the excited state couples to phonons, 
but our results could be easily generalized if phonon cou- 
pling was also present in the dot groundstate^ 

We next employ a density matrix approach .^'^^'^^'^^'^ 
To this end, we start from the density matrix pij = {(Jij)^ 
where the brackets denote tr(pcrij), and compute its 
equation of motion from the Heisenberg equation of mo- 
tion 



^eff^^j) + {[(^ij^(^ee]Bee) • (28) 



Because of the phonon coupling, the density matrix cou- 
ples to the phonon- assisted density matrix^'^ (aijbk)- If 



TABLE I: Correlation functions used in our density matrix 
approach. We consider all correlation functions which contain 
less than three phonon operators. The equations of motion for 
the different correlation functions are given in Appendix [A] 

Correlation function Symbol Expression 



Density matrix 
Phonon amplitude 
Phonon-assisted density matrix 
Two-phonon amplitude 
Phonon occupation 
Two-phonon-assisted density matrix 



pij 

Sk 

Skk' 
rikk' 
Sij^kk' 
nij^kk' 



ih) 

{{Kh')) 

{{^^JK^k') 



we compute next the equation of motion for the phonon 
assisted density matrix, we find that it couples to a 
density matrix with two phonon operators, whose time 
evolution, in turn, involves density matrices with three 
phonon operators. This refiects the fact that the inter- 
acting dot-phonon system constitutes a highly non-trivial 
many-body problem, with an infinite hierarchy of equa- 
tions of motion for the density matrices. To truncate 
this infinite hie rarchy, it is convenient to employ a cumu- 
lant expansion^^^^^ where one introduces correlation 
functions by subtracting from the higher-order density 
matrices contributions of lower-order density matrices. 
For instance, we obtain for the phonon assisted density 
matrix 



h)) , 



(29) 



where ((•)) denotes the correlation function. On general 
physical grounds, we expect that higher-order phonon- 
assisted correlation functions will play no significant role 
in the dynamics of pij. 

In our density matrix equation approach we first select 
a set of representative correlation functions. Herewe re- 
strict ourselves, in accordance to related work,SSliII to 
those correlation functions which contain less than three 
phonon operators (see Table [l] for a complete list). Their 
equations of motions are obtained by using the Heisen- 
berg equation of motion, performing the cumulant expan- 
sion, and keeping only the relevant correlation functions. 
By neglecting higher-order correlations, we truncate the 
infinite hierarchy of equations of motions, and obtain a 
closed set of equations of motion (see Appendix [A|) . This 
set can be solved through direct numerical integration. 
Results of these calculations, and a comparison with the 
perturbative approach, will be given in the next section. 



III. RESULTS 
A. Model 

In our simulations we use a quantum dot and phonon 
model similar to Ref. ^20^, with Gaussian electron and hole 



TABLE II: Material, dot, and cavity parameters used in our 
simulations. The phonon parameters are representative for 
GaAs. For the electron and hole wavef unctions, we assume 
Gaussians with a full width of half maximum (FWHM) of Liat 
and Lz along the lateral and growth direction, respectively. 
Throughout we use an excit on- cavity coupling strength of 150 
jieV and a radiative dot lifetime of 7 ns. The cavity quality 
factor is Q = 15000 unless stated differently. 



Parameter 



Symbol Value 



mass density p 5.37 g/cm~^ 

sound velocity q 5110m/s 

deformation potential for electrons De —14.6 eV 

deformation potential for holes Dh —4.8 eV 

in-plane confinement (FWHM) Liat 10 nm 

confinement in 2;-direction (FWHM) Lz 4 nm 

biexciton binding energy Ab 2 meV 

radiative decay time of dot Trad 7 ns 

exciton-cavity coupling g 150 /xeV 

cavity photon energy cjcav 1-3 eV 
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FIG. 2: (Color online) Cavity feeding rate for the phonon- 
assisted transition from the exciton to the cavity. We use 
the dot and material parameters listed in Table. [H] and an 
exciton-cavity coupling g = 150 /ieV. The solid and dashed 
lines report results obtained from Eqs. (20) and (13) for the 



independent Boson model (IB) and the Schrieffer-Wolff (SW) 
perturbation approach, respectively. 



wavefunctions, and phonon parameters representative for 
GaAs (see Table |ll). Quite generally, we expect that a 
more realistic description scheme, including effects due 
to strain, the more complicated valence band structure, 
or the heterogeneous material composition, might some- 
what alter our results. On the other hand, our model has 
proven successful i n repro ducing the dominant effects of 
phonon dephasing,^MESI ^nd is sufficiently generic to 
include possible modifications from a microscopic model 
description by simply adapting the effective parameters 
(quantum dot form factor, phonon coupling). We only 
consider deformation potential phonon coupling, which 
is expected to be the dominant contribution. For the 
cavity feeding through decay of the exciton, the exciton 
phonon coupling is of the fornPSEH 



2pci 



ikr 



(30) 

where k denotes the wavevector, and ^e,/i(^) are the elec- 
tron and hole wavefunctions for which we assume Gaus- 
sians (see Table |ll| for parameters). Piezoelectric cou- 
pling, described by performing angular aver ages o ver the 
longitudinal and transverse modes separatelyj ^^*^^ * turned 
out to give only negligible corrections, and was thus ne- 
glected. Coupling to LO phonons leads to the formation 
of a tightly bound polaron.l^^I^The resulting energy shift 
can be easily absorbed into the exciton and biexciton en- 
ergies. The polaron-mediated relaxation channel for ex- 
cited states, provided by the anharmonic decay of LO 
phonons, is expected to be of minor importance since we 
are dealing with the exciton and biexciton ground states 
(where polarons are stable) and are only interested in en- 
ergy transfers much smaller than the LO phonon energy. 



B. Phonon-assisted cavity feeding 

Figure [2] shows the cavity feeding rates for the phonon- 
assisted transition from the exciton to the cavity (and the 
dot ends up in the groundstate). The solid lines report 
results from the independent Boson model, Eq. ([20]), and 
the dashed lines results of Eq. ([l3| for the Schrieffer- 
Wolff perturbation approach. At low temperatures, the 
feeding rate T^g as a function of detuning is strongly 
asymmetric. For positive detuning A = — Eg — c^cav, 
i.e., when the exciton has a higher energy than the cav- 
ity mode, the transition from the exciton to the cavity 
is accompanied by a phonon emission. In this regime 
there is a substantial scattering probability even at low 
temperatures. With increasing temperature the phonons 
become thermally populated, as described by the Bose- 
Einstein factors fik in the scattering rates of Eqs. (20) 



and Eq. (13), leading to a noticeable scattering proba- 
bility also for negative detunings. These processes are 
associated with phonon absorption. At sufficiently high 
temperatures, say above 20 K, the feeding rate as a func- 
tion of detuning becomes almost symmetric. 

The dependence of the feeding rate on A is governed 
by the exciton form factor, as given by the integral in 
Eq. (30). For the Gaussian carrier wavefunctions, the 



width of the Fourier-transformed probability distribution 
is again a Gaussian, with a width given by the inverse of 
the confinement length L. Due to energy conservation 
in the scattering, the wavenumber /c, at which the form- 
factor has to be computed, is determined by cik = A. 
Thus, a small confinement length L translates to a broad 
/c-space distribution of the formfactor, and, in turn, to a 
broad re^(A) distribution. In contrast, for a larger con- 
finement length the feeding rate as a function of detuning 
becomes narrower. Additionally, the feeding rate is di- 
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FIG. 3: (Color online) Spectra for driven dot-cavity system 
and for different quality factors. We use a level scheme con- 
sisting of the quantum dot ground and excited states, and 
the cavity. The panels report results for a detuning of (a) 
A = -50 /ieV, (b) A = 0, and (c) A = 50 /ieV. The temper- 
ature is set to T = 10 K. 



rectly proportional to the Purcell enhancement^^ (^/A)^ 
caused by the photon confinement within the microcav- 
ity. We emphasize that the scattering rates of both the 
independent Boson model and the Schrieffer- Wolff per- 
turbation approach directly scale with and do not 
depend on the quality factor Q of the cavity. However, 
Q determines the lifetime of the cavity photon, and thus 
has a decisive impact on the efficiency of phonon-assisted 
cavity feeding. 

From the comparison of the solid and dashed lines in 
Fig. [2] we observe that the results of the Schrieffer- Wolff 
perturbation approach and the independent Boson model 
are in almost perfect agreement, for temperatures below 
say 20 K. For elevated temperatures, the Schrieffer- Wolff 
approach overestimates the feeding rate. We attribute 
this to the neglect of re-absorption processes of phonons, 
which are naturally included within the independent Bo- 
son model. From now on we will only consider the feed- 
ing rates for the independent Boson model, Eq. (20), al- 



though our results would look very similar for T^g com- 
puted within the Schrieffer- Wolff perturbation approach. 



C. Spontaneous emission lifetime 

We next consider the situation where the dot is initially 
in the exciton state, and subsequently decays through ei- 
ther coupling to the leaky cavity mode, phonon-assisted 
cavity feeding, or radiative decay. From Fig. [3| which 
shows luminescence spectra for a driven dot-cavity sys- 
tem, to be discussed further below, we observe that for 
the chosen dot and cavity parameters the system operates 
in the strong coupling regime for quality factors above 
approximately Q = 2000. Indeed, for Q — 5000 the de- 
cay transients in Fig. [4] display at early times population 
oscillations between exciton and cavity (best visible for 
zero or small detunings). These oscillations are a clear 
signature of strong coupling. 

The solid lines report results of our density matrix ap- 
proach. When comparing the decay for positive and neg- 
ative detunings ±A, we find that the decay is always 
faster for positive detunings than for negative detunings. 




FIG. 4: (Color online) Exciton decay for a dot embedded 
in the cavity, and in presence of phonon-assisted transitions. 
The solid and dashed lines report results of our density ma- 
trix and master equation approach, respectively. For clarity, 
the decay transients for different detunings are offset in time. 
Throughout the decay for positive (blue) detuning (exciton 
energy larger than cavity energy cJcav) is faster than for neg- 
ative (red) detuning. In all simulations we use Q = 5000 and 
T = 10 K. The inset reports results for the correlation func- 
tion C{t) of the independent Boson model, Eq. ( 19 ). The solid 
lines correspond to the analytic expressions, and the symbols 
show the results of the density matrix approach. Both results 
are in perfect agreement, thus demonstrating the accuracy of 
our numerical approach. 



This is due to the enhancement of phonon emissions in 
comparison to phonon absorptions, attributed to the dif- 
ferent bosonic factors of n/e + 1 and flk in the scattering 
rates. The dashed lines show results of the correspond- 
ing master equation simulations. We find nice agree- 
ment, with exception of small detunings, where probably 
higher-order scattering or renormalization processes play 
some role. 

In the decay of the exciton, the population decays 
mono-exponentially at times above approximately 20 ps. 
We can thus assign a lifetime to the excited quantum dot 
state, which we shall denote as the spontaneous emission 
(SE) lifetime. In Fig. [s] we show the extracted lifetimes 
rsE{^) for a variety of quality factors and temperatures. 
Quite generally, we find that rssi^) dramatically de- 
creases when exciton and cavity come into resonance. 
This is attributed to the formation of a polariton with 
a significant cavity admixture, such that the exciton can 
decay through the leaky cavity, rssi^) has a broad dip 
for small quality factors, which significanty narrows for 
higher Q values. 

The gray lines in Fig. |5] show results of simulations 
where phonon scatterings have been artificially neglected. 
When phonon scatterings are included, tse as a func- 
tion of detuning becomes strongly asymmetric. In par- 
ticular for the higher quality factors, rsE{^) is signif- 
icantly reduced at positive detunings. This is due to 
the aforementioned cavity feeding through phonon emis- 
sion, which is the dominant scattering mechanism at low 
temperatures. With increasing temperature tse{^) be- 
comes also reduced at negative detuning s, as a result of 
thermally activated phonon absorption. In the insets of 
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FIG. 5: (Color online) Spontaneous emission (SE) decay time of initially excited quantum dot for different quality factors and 
temperatures (same line labeling as in Fig. [2]). The gray lines report results where phonon scatterings have been artificially 
neglected. In the insets we compare, for a smaller range of detunings, the density matrix approach (symbols) with the 
perturbation approach (solid lines), finding perfect agreement throughout. 



the figure we show, for a smaller range of detunings, re- 
sults obtained within our density matrix approach (sym- 
bols) and master equation approach (solid lines). One 
clearly observes that both approaches give almost iden- 
tical results, for all the quality factors and temperatures 
investigated. We thus conclude that the master equation 
approach with phonon feeding rates, computed from ei- 
ther Eq. (p^Sl) or (l20|), provides an accurate description 



scheme for the problem under consideration. From now 
on we will only use the master equation approach, de- 
scribed in Sec. |IID[ which provides a simple and versatile 
simulation tool. 



branch, which is obtained by integrating in the spectra 
over the spectral region indicated by the dashed lines. 

We observe, in accordance to related work,^^ a highly 
non-trivial dependence of the polariton intensity as a 
function of exciton-cavity detuning. To understand this 
dependence, we first note that the system is driven 
through the quantum dot state, and the transition with 
the higher intensity is the one associated with the quan- 
tum dot state, unless there is an efficient feeding of the 
cavity. Quite generally, the polariton modes can be ap- 
proximately determined from the eigenvalues and eigen- 
vectors of the effective two-mode hamiltoniari^ 



^eff 



D. Spectrum 

We next discuss the influence of phonon-assisted cavity 
feeding on the spectral properties of the coupled dot- 
cavity system. In our simulations we consider the level 
scheme depicted in Fig.jlja), consisting of the dot ground 
and excited state, and the cavity mode. The system is 
weakly driven by an external laser pulse, as described by 
a Lindblad operator which promotes the system from g to 
e. If the pump rate is sufficiently small, in our case much 
smaller than the radiative decay rate of the quantum 
dot, only the overall intensities of the spectrum (but not 
the relative intensities of the upper and lower polariton 
branches) depend on the pump rate. In our simulations 
we also include a dephasing rate for the exciton, primarily 
for visualization purposes to broaden the exciton line. 
We found that for detunings larger than the polariton 
linewidths this dephasing has an only minor impact on 
our results. 

Figure [oj^b) shows the computed spectra for different 
exciton-cavity detunings. At zero detuning we observe 
an anticrossing of the lower and upper polariton branch, 
which is the signature of strong coupling in the frequency 
domain. In panel (a) of the figure we plot the 
intensity of the lower (red) and upper (blue) polariton 



(31) 



where the diagonal terms account for the cavity and ex- 
cited dot states, including radiative damping, and the 
off-diagonal terms for the dot-cavity coupling. In the 
following we discuss for the plait pattern of the polari- 
ton intensities, depicted in Fig.|6ja), the different regions 
indicated by numbers in circles. (1) For large negative 
detunings, the lower polariton branch has a dominant 
exciton character, and the probability for cavity feeding 
through phonon absorption (which is needed since the 
exciton has a lower energy than the cavity) is negligible. 
Here, the exciton decay becomes enhanced through the 
cavity, but the system decays at the exciton frequency. 
In this regime the lower polariton branch has the higher 
intensity. (2) When the detuning is reduced, phonon- 
assisted feeding becomes more likely, due to the small 
energy mismatch, and the upper polariton (which has a 
predominant cavity contribution) emits with a slightly 
higher intensity. We will show in a moment that this 
cavity transfer through phonon absorption can be sup- 
pressed at even lower temperatures. (3) When the de- 
tuning is further reduced, the lower polariton branch ac- 
quires a larger contribution from the cavity, and can thus 
decay more efficiently through the leaky cavity. In this 
regime the lower polariton has a higher intensity. (4) At 
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FIG. 6: (Color online) Spectrum of driven dot-cavity system, 
as depicted in Fig. [l] (b) Spectrum as a function of detuning 
between exciton and cavity. One observes the formation of 
polariton states, which anticross at zero detuning as a result 
of strong coupling. In our simulations we use a quality factor 
of Q = 15000, a temperature of T = 10 K, a radiative dot 
lifetime of Trad = 7 ns, and an exciton dephasing time of 
50 ps. (a) Intensity of lower and upper polariton branch. 
The intensity is integrated over the regions indicated by the 
dashed lines in panel (b). The bright lines show results of 
simulations where phonon scatterings have been artificially 
neglected. The numbers in circles are used for the discussion 
in the text, and the dashed line indicates the position of zero 
detuning where the two polariton modes have equal intensity. 



zero detuning the two polariton modes have equal exci- 
ton and cavity character, and the intensities of the two 
polariton branches are equal. (5) For positive detuning, 
the upper (blue) polariton branch has a stronger exci- 
ton character. When the detuning is sufficiently small, 
the upper polariton has a noticeable cavity contribution, 
and thus emits with higher intensity (analogously to the 
case of red detuning). (6) At larger positive detunings 
phonon-assisted cavity feeding through phonon emission 
becomes the dominant mechanism, and the lower (cavity- 
like) polariton has a much higher intensity. 

From the figure we clearly observe a strong asymme- 
try between negative and positive detunings, which is 
attributed to the different probabilities for cavity feed- 
ing through phonon absorption or emission. These re- 
sults are in complete accordance to the above discussion 
about the spontaneous emission lifetime. Simulations 
performed by artificially neglecting phonon scatterings 
(bright lines) only give a Purcell enhancement around 
zero detuning. Otherwise the results totally differ from 
those of the simulations including phonon-assisted cavity 
feeding, thus highlighting the importance of this scatter- 
ing channel. 

In Fig. [T] we investigate the polariton intensities for 
different quality factors and temperatures. Panel (a) 
reports the influence of the quantum- dot- like polariton 
(lower polariton branch for negative detunings, and up- 
per branch for positive detunings) intensity on the quality 
factor (we use T = 10 K). For small Q values, say be- 
low 2500, this branch is strongest for all detunings. For 
larger Q values the effect of cavity feeding becomes more 
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FIG. 7: (Color online) Relative intensity of quantum- dot- like 
polariton branch as a function of (a) quality factor Q (with 
T = 10 K) and (b) temperature (with Q = 15000). The 
dashed lines indicate zero detuning. The relative intensity 
of the polariton modes is computed by integrating over the 
area for the lower (red) [upper (blue)] polariton for negative 
[positive] detunings, as indicated in Fig. |6jb), and dividing 
by the sum of the upper and lower polariton intensity. For 
discussion see text. 



important, and the intensity significantly drops for posi- 
tive detunings where the exciton can decay into a cavity 
photon. We observe that the plait pattern intensity de- 
pendence as a function of detuning is most pronounced 
for high Q values. Here, the cavity photon has the longest 
lifetime, and accordingly phonon-assisted cavity feeding 
becomes most efficient. Nevertheless, the influence of Q 
on the polariton intensity is not overly pronounced in this 
strong-coupling regime. Panel (b) reports the relative 
intensity of the lower polariton branch as a function of 
temperature. The values at 10 K exactly corresponds to 
the graph of Fig. [6j which we have discussed above. We 
observe that cavity feeding through phonon absorption 
(point 2 of the above discussion) becomes completely sup- 
pressed at low temperatures. On the other hand, feeding 
gains importance at elevated temperatures. At tempera- 
ture above 30 K the system already predominantly emits 
at the cavity frequency. 

We note that the plait pattern of the intensity also 
depends decisively on other feeding channels of the cav- 
ity. Such feeding is usually present in experiments due 
to cavity-assisted relaxatio n pro cesses of carriers initially 
excited in the wetting layer We model it by an addi- 
tional Lindblad operator that brings the system directly 
from g to the cavity (see Fig.jl]). Our results (not shown) 
reveal that the characteristic plait pattern almost disap- 
pears when the cavity feeding rate is twice as large as the 
dot pumping rate. Thus, the polariton intensity depen- 
dence could serve as a sensitive measure for additional 
cavity feeding processes. 



E. Biexciton decay and driven biexciton system 

In the following we consider a more complete level 
scheme for the quantum dot, depicted in Fig. [TJb), con- 
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FIG. 8: (Color online) Time evolution of biexciton and exci- 
ton population in the biexciton cascade decay. In the simu- 
lations, at time zero only the biexciton is populated, which 
decays via the exciton states to the ground state. The dif- 
ferent rows report the population transients for the exciton 
Xfast , that couples with g cos strongly to the cavity, the ex- 
citon Xsiow, that couples with ^sin^ weakly to the cavity, 
and the biexciton XX. The different columns report results 
for different mixing angles 6 [see Fig.[T]^b)]. In our simulations 
we use Q — 15000 and T = 10 K. The contour lines indicate 
where the population has decayed to 37% of its maximum 
value. For the fast exciton and the biexciton decay these con- 
tour lines are also magnified to make the features at early 
times better visible. The dashed lines indicate the position 
where the biexciton is in resonance with two cavity photons. 



sisting of the dot groundstate, two exciton states x and 
y with orthogonal, linear polarization, and the biexciton 
state XX. We also include states where, in addition, one 
or two photons are present in the cavity. We will pri- 
marily investigate how our conclusions for the two-level 
system, discussed above, become modified when consid- 
ering this more complete dot description. 

From the beginning we introduce a few simplifications. 
First, we neither consider (spin forbidden) dark states, 
which might play some role for small detunings or for 
the driven biexciton system, nor charged exciton states, 
which seem to be ubi quitous for quantum dots embedded 
in nanocavities.'^^'^ Also the finestructure splitting be- 
tween the two exciton states is ignored. In principle, such 
splitting could be easily included in our simulations, but 
we do not expect that the resulting modifications are of 
relevance for the relatively large detunings (up to several 
meV) considered here. In our simulations we assume that 
both exciton states x and y can couple to the cavity, and 
introduce a mixing angle 9 that determines the relative 
coupling strengths of gcosO and ^sin^, respectively. We 
will only be interested in small mixing angles, and will 
consequently denote the excitons with the fast {gcosO) 
and slow {g sin 0) decay characteristics as Xfast and Xgiow, 
respectively. 
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FIG. 9: (Color online) Steady state exciton and biexciton 
populations of a driven biexciton system. We use Q = 15000 
and T = 10 K. The labeling of the exciton and biexciton 
states is identical to Fig. [S] The different columns correspond 
to different pumping rates, indicated on top of the panels. 



Figure [S] shows simulation results of a dot that is 
initially populated with a biexciton, and subsequently 
decays through radiative couplings and phonon-assisted 
transitions to the cavity. No pumping is considered here. 
The situation under investigation approximately corre- 
sponds to experiments where at time zero electron-hole 
pairs are excited in the wetting layer, or in excited dot 
states, and consecutively relax on a short timescale to 
the biexciton state of lowest energy. Let us first consider 
a zero mixing angle (left column). In accordance to the 
spontaneous emission lifetimes of the exciton, shown in 
Fig. [sj we observe for the biexciton (lower left panel) 
an asymmetric detuning dependence for the population 
decay, which, however, is less pronounced than for the 
exciton. The fastest decay is now at a detuning corre- 
sponding to the biexciton binding energy A5 of 2 meV, 
where the biexciton to exciton decay is in resonance with 
the cavity. 

For our quantum dot model, the phonon-assisted cav- 
ity feeding rates for the biexciton are the same as for the 
exciton. According to the discussion of Sec. |n| we have 
to use in the scattering rates the difference Xee,k — ^gg,k 
of biexciton and exciton phonon coupling constants. If 
we neglect possible modifications of the carrier wavefunc- 
tions in the biexciton state, which is certainly a good 
approximation that will be adopted in this work, the 
scattering rates for phonon-assisted biexciton and exci- 
ton decay are the same, apart from the inclusion of the 
biexciton binding in the detunings. 

A striking feature of the biexciton decay depicted in 
Fig. |8] is the sharp feature around 1 meV, best visible 
from the contour line that indicates where the popula- 
tion has dropped to 37% of its initial value (see closeup 
where the time scale has been reduced by a factor of 50) . 
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The detuning of 1 meV corresponds to the two photon 
resonance condition 



(32) 



where the biexciton ca n dec ay by emitting two cavity 
photons simultaneouslyJ^^'^ Here is the exciton en- 
ergy. Because the transition involves a tunnehng-hke 
transition, through the off-resonant exciton state, the de- 
tuning width where this process is efficient is extremely 
narrow. 

Through the biexciton decay, the exciton states be- 
come populated and consecutively further decay to the 
groundstate. Let us first concentrate on the fast exciton 
component Xfast- We observe for the population decay a 
strongly asymmetric detuning dependence, in accordance 
to our previous discussion for the two-level dot model. 
However, the minimum is now shifted to larger detunings. 
This is due to the fact that the biexciton decays fast, and 
correspondingly the exciton population builds of quickly. 
Thus, the exciton decay reflects the combined effects of 
population buildup, through the biexciton decay, and ex- 
citon decay. Incidentally, the fastest exciton decay is at 
the detuning where the biexciton is in resonance with two 
cavity photons, and also the exciton decay is fast as a re- 
sult of the Purcell enhancement. We note that for this 
detuning also the population that is channeled through 
the exciton is smallest, since the biexciton can bypass the 
excitons in the direct two-photon decay. 

We next investigate the decay of the slow exciton com- 
ponent Xsiow For zero mixing angle, ^ = 0, exciton 
feeding (through the biexciton) and decay are not en- 
hanced by the cavity. We observe in Fig. |8] that for suf- 
ficiently large negative detunings, where the biexciton 
decays without a significant Purcell enhancement, both 
fast and slow exciton components are populated equally. 

decays with the radiative decay time Trad, in con- 
trast to the fast component, whose decay is enhanced 
through cavity coupling. For increased detunings, the 
biexciton decay becomes enhanced through phonon feed- 
ing and the Purcell effect, and decays primarily via the 
Xfast channel. Thus, for ^ = the decay of the slow 
exciton component directly mirrors the decay character- 
istics of the biexciton. In presence of a non-zero mixing 
angle Xgiow acquires a small cavity admixture, through 
the coupling ^sin^, and its feeding and decay become 
enhanced. Indeed, we observe a pronounced lifetime de- 
crease for small cavity detunings, with the minimum be- 
ing at the two-photon resonance. 

Finally, in Fig. [9] we analyze the steady state popula- 
tions of the driven biexciton system depicted in Fig. [l] 
We chose pumping rates, indicated on top of the panels, 
which are smaller or larger than the radiative decay rates 
of the isolated exciton states. For essentially all pumping 
rates we observe a significant population of the biexciton 
state. This is due to the fact that the Xgiow population 
is usually large, due to the long lifetime originating from 
the weak cavity coupling. Thus, there exists an efficient 
population channel of the biexciton via the intermediate 
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FIG. 10: (Color online) Simulated spectra for driven biexciton 
system. We use Q — 15000, T = 10 K, a mixing angle of 
^ = 5°, and a pump rate of 1/(5 ns). To make the biexcitonic 
features better visible, we use a logarithmic color scale. The 
intensities of the polariton modes are computed in accordance 
to the prescription given in the caption of Fig.|6] The numbers 
in circles are used in the discussion in the text. 



Xgiow state. The steady state populations of all states 
reflect the respective decay characteristics. The popula- 
tions are large when the lifetime is long, and small when 
the lifetime is short. As expected, with increasing pump- 
ing rate the steady state populations increase. 



F. Biexciton spectra 



We finally compute luminescence spectra for a driven 



biexciton system. Figure [TO] shows the spectra as a func- 
tion of detuning. We use a logarithmic color scale to 
make the weak biexciton features better visible. In the 
luminescence spectra we observe, in addition to the lower 
and upper polariton branches of the coupled exciton- 
cavity system, peaks that are shifted by approximately 2 
meV to the red, as indicated in the figure by (1). These 
features are associated with the biexciton to exciton de- 
cay. (2) The anticrossing at zero detuning is attributed 
to the biexciton decay to the strongly coupled upper 
and lower exciton-polaritons.^^ At point (3) we observe 
a crossing with the lower exciton-polariton branch. As 
the two branches are associated with the cavity-like po- 
laritons, this crossing point precisely indicates the point 
where the biexciton is in resonance with the two-photon 
transition. (4) Finally, when the biexciton to exciton 
transition is in resonance with the cavity mode, we ob- 
serve an anticrossing of the polaritons associated with 
this transition. At the same spectral position there is 
additonally the emission of the cavity-like polariton of 
the coupled exciton-cavity system. 

The intensity of the lower (red) polariton branch, 
shown in panel (a), exhibits a small bump at the point 
of two-photon resonance. When comparing the results of 
Fig. [lOja) with those of the two-level system. Fig. [6]^ a), 
we find that the relative intensities are very similar, and 
only the detuning dependence of the absolute intensities 
differ. This is due to the more complicated excitation 
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dynamics in case of biexciton pumping. Nevertheless, 
the main conclusions for the dot-cavity system regarding 
quality factor and temperature prevail. 



IV. SUMMARY AND CONCLUSIONS 

To summarize, we have investigated phonon-assisted 
transitions for a coupled quantum dot-cavity system, 
where the exciton decays into the cavity and the energy 
mismatch is compensated through phonon emission or 
absorption. Such processes play an important role for 
energy detunings of the order of a few meV, and for dot- 
cavity systems operating in the strong coupling regime. 
The basic mechanism underlying phonon-assisted cavity 
feeding can be encapsulated in the effective Hamiltonian 



photon sources, '^n^E^'I^ or in quantum information sci- 



H' 



3 ( 



(33) 



which has been derived within the framework of a 
Schrieffer- Wolff transformation [Eq. (pT|]. It describes a 
combined phonon and cavity effect, where Ogf^a^ accounts 
for the exciton decay and the creation of a cavity photon, 
and the phonon field operators 6^ and for the accompa- 
nying emission or absorption of a phonon. Although such 
processes can in principle also occur for quantum dots 
outside a microcavity, they become strongly enhanced 
within a cavity, because of the Purcell enhancement of 
<7/A. 

In addition to the Schrieffer- Wolff perturbation ap- 
proach, we have derived a similar expression for the inde- 
pendent Boson model [Eq. (20)], where the phonon dy- 
namics has been accounted for in an exact manner. We 
have also analyzed the problem within a numerical den- 
sity matrix approach, without employing any perturba- 
tion procedure. For realistic dot and cavity parameters, 
all approaches give more or less the same results, at least 
for temperatures below 20 K. 

We have investigated the role of phonon-assisted cav- 
ity feeding in the exciton and biexciton decay, and have 
found a strong detuning asymmetry in agreement with 
experiment. This has been attributed to the different 
importance of phonon emissions and absorptions. For the 
spectra of a driven dot system, we have found an intrigu- 
ing plait pattern for the intensities of the upper and lower 
polariton branches. We have shown that these features 
prevail for more realistic dot level schemes. For the biex- 
citon cascade decay, we have analyzed the different decay 
channels via the intermediate exciton states, and have 
found a fast and slow decay component, in agreement 
with experiment. For a driven biexciton state, an effi- 
cient population channel of the biexciton via one of the 
exciton states (weakly coupled to the cavity) has been 
identified. 

Phonon-assisted cavity feeding plays a dominant role 
in cavity- QED experiments whenever the dot and cavity 
are detuned by only a few meV, and might be of impor- 
tance for quantum dot-based lasers,!^ single or entangled 



ence." 
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APPENDIX A: DENSITY MATRIX APPROACH 

In this appendix we give some details of our density 
matrix approach. For each correlation function ((A)) we 
derive the dynamic equation from the Heisenberg equa- 
tion of motion 

^|(^) = (^J^iff - ^^eff^) + ([A^eelSee) , (Al) 

and subtract the contributions from the lower-order den- 
sity matrices. In the following we only consider the non- 
trivial terms according to phonon coupling. The contri- 
butions due to i^eff can be added in a straightforward 
fashion. 

As for the density matrix of the two-level system we 
obtain 



^Pge — ^Pge ^ge 
'^Peg — ^Peg ^eg •> 



(A2a) 
(A2b) 



together with pgg = pee =0. Here we have introduced 
the abbreviations 



E 

k 

E 

k 



Xee,k{sk^sl) (A3a) 

>^ee,k{Seg,k^sle,k) (A3b) 



for the energy renormalization as well as the source 
terms Qeg and = 1^*^, which describe the coupling 
to the phonon-assisted density matrices. For the time 
evolution of the phonon density matrices we get 

iSk = Xee,kPee (A4a) 
iSkk' = Xee,kSee,k' Xee,k'See,k (A4b) 
i'Tlkk' = Ke,kKe,k' ~ '^ee,k'See,k J (A4c) 

and for that of the phonon-assisted density matrix 

^^ge,k — ^^ge,k H~ '^ee,/c(l Pee)Pge H~ ^kPge H~ ^ge,k 
'^^eg^k ~ ^^eg,k Xee,k Pee Peg ^kPeg ^eg,k 

'^^ee^k — ^ee,fe(l Pee)Pee 



'^^gg^k 



0. 



(A5) 



We have introduced for the coupling to the higher-order 
phonon correlations the abbreviations 



^k = Xee,k' {Skk' + ^kk 

k' 

^eg,k — ^ ^ Xee,k' (^^eg,kk' H~ ^eg, 



kk' 



(A6a) 
(A6b) 
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For the two-phonon- assisted density matrices we finally together with Sgg^kk' = '^gg,kk' ~ ^' "^^^^ provides us with 
obtain a closed set of equations. In our computational approach, 

we introduce a surrogate hamiltonian, with a finite num- 
k'Sge.k ber of representative phonon modes (typically a few hun- 

+ ^ee,/e ((1 " Pee)Sge,k' - PgeSee.k') dred),'^^ and solvc the cquatious of motion through direct 

. ' n-^ \ \ / A rr \ numerical integration. 

+ Aee,/c' ((1 - PeeJSge^k " Pge^ee^k ) (A7a) 
^^eg,kk' — ^^eg ,kk' ^k^eg ,k' ^k' ^eg,k 
^ee,k {^Pee^eg ^k' Peg^ee^k'^ 

- Aee,fe' {peeSeg^k + PegS^e.k ) (A7b) 
'i^See.kk' = ~ '^Pee){^ee,kSee,k' + ^ee,k'See,k) (A7c) 

and 

'^'f^ge^kk' ~ ^"^ge^kk' ^k^eg,k' ~^ ^k' ^ge,k 

+ ^ee,k ((1 ~ Pee)slg,k' ~ PgeKe,k') 

+ Ke,k' ( Pee Sg^^j, + PgeSee,k ) (A8a) 

'^'^eg,kk' — ^^eg,kk' ^k^ ge^k' ^k' ^^g ,k 

~ ^ee,k ( Pee ^le,k' ~^ PegKe,k') 

- Ke,k' ((1 - Pee)Seg,k " Peg^ee^k ) (A8b) 
'i^ee^kk' = (1 " 2/)ee) (Aee,fc5ee,/c' ~ Aee,/c' 5ee,/c) , (A8c) 
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